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ABSTRACT 

The general formalism to describe single spin asymmetries in hadron-hadron 
high energy and large pt inclusive production within the QCD-improved parton 
model and assuming the factorization theorem to hold for higher twist contri- 
butions is discussed. Non zero values of the single spin asymmetries originate 
from and reveal non perturbative universal properties of quarks: the quark dis- 
tribution and fragmentation analysing powers, which need not be zero, provided 
the quark intrinsic motion is taken into account. A simple model is constructed 
which reproduces the main features of the data on the single spin asymmetries 
observed in inclusive pion production in pp collisions. 



According to the QCD factorization theorem the differential cross-section for the 
hard scattering of a polarized hadron A with spin Sa off another polarized hadron 
B with spin Sb-, resulting in the inclusive production of a hadron C with energy Eq 
and three-momentum p^, A,Sa + B,Sb^C + X, can be written as 0Si 



X pl^^i^ U,a{x.) p^/^f f„s{x,) M.... Mv,A.;A' A' ^A^l'll^) 



* Talk delivered by M. Anselmino at the Adriatico Research Conference on "Trends in Collider Spin 
Physics", Trieste, Italy, December 1995, and at the International Conference on Hadron Structure: 
High-Energy Interactions: Theory and Experiment (Hadron Structure 96), Stara Lesna, Slovak 
Republik, February 1996. 
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where fi/i{xi) is the number density of partons i with momentum fraction Xi inside 
hadron I {i = a,b; I = A,B) and p^l^'y'{xi) is the hehcity density matrix of parton 

i ' i 

i inside the polarized hadron I. The Ma^,a^;A„,Aj^'s are the hehcity amphtudes for the 



elementary process ab cd; if one wishes to consider higher order (in as) contribu- 

A A' 

tions also elementary processes involving more partons should be included. D^"' y (z) 
is the product of fragmentation amplitudes 

^'cJo=^,,,,^^x.Ac:A.^^I.,A',;A,, (2) 

where the T stays for a spin sum and phase space integration of the undetected 

JX,Xx 

particles, considered as a system X. The usual unpolarized fragmentation function 
Dc/c{z), i-S- the density number of hadrons C resulting from the fragmentation of 
an unpolarized parton c and carrying a fraction z of its momentum is given by 

Dc/c{z) = \ E Dl'^X^z). (3) 

For simplicity of notations we have not shown in Eq. (|l]) the scale dependences 
in / and D] the variable z is related to Xa and x^ by the usual imposition of energy 
momentum conservation in the elementary 2 — > 2 process §, z = —{xat + Xbu)/xaXbS 
where s, t, u are the Mandelstam variables for the overall process AB — »• CX whereas 
s, t, u are those for the elementary process ah — > cd. The corresponding amplitudes 
are normalized so that 

da 1 1 ^ I " 12 , X 

d£ = 16^4, • (4) 

Eq. (|T|) holds at leading twist and large pt values of the produced pion; the 
intrinsic fc^ of the partons have been integrated over and coUinear configurations 
dominate both the distribution functions and the fragmentation processes; one can 
then see that, in this case, there cannot be any single spin asymmetry. 

Suppose, in fact, that hadron B is not polarized so that Eq. (|l]) reads 

d'Pc ~ .LA.,A,;A,tA,A,;A,2i VT. IQ.P 

X Ptt /'^m(^«) /Vb(^'>) ^A.,A,;A.,A, <,A,;A'.,A, DI^^z) . 

Then total angular momentum conservation in the (forward) fragmentation pro- 
cess [see Eq. (H)] implies = A^; this, in turns, together with helicity conservation 
in the elementary processes, implies A^ = A^. If we further notice that, by parity 
invariance, X^a -^a'" a'' — ^cic independently of Ac and that 

Y, |iVW,vP = 2(W);|, (6) 



dt 
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we remain with J2\ P\^x^ = 1- Moreover, in the absence of intrinsic k± and initial 
state interactions, the parton density numbers fa/Ai^a) cannot depend on the hadron 
spin and any spin dependence disappears from Eq. (|1|), so that 

2^ J Ia/A[Xa) Jb/B[Xb) 



C a,b,c,d 



This imphes that the single spin asymmetry, 

^^A^ B^CX^p^^ _ ^^A^B^CX^_p^^ 



A 



N 



da^'B^^^iPj) + da^'B-^^^i-pj) ' 

is zero, t (i) refers to hadron A spin orientation up (down) with respect to the 
production plane; any other spin orientation with no up or down component would 
result in a vanishing asymmetry due to parity invariance. 

Eq. (|l]) can be generalized with the inclusion of intrinsic k± i and this can avoid 
the above conclusion; for example i, the observation of a non zero fc^ of a final 
particle C with respect to the axis of the jet generated by parton c does not imply 
any more = and allows a non zero value of the asymmetry through the quark 
fragmentation analysing power 

Dc/q,s^ {Z, ki_) + Dc/q,-s, {Z, fcx) 
Dc/g,s,{z, fcx) - Dc/g,s,{z, -fc^) 
Dc/q,s,{z, fcx) + Dc/q,s,{z, ■ 

A non zero value of the above quantity - the analysing power or single spin 
asymmetry for the quark fragmentation process - is allowed by parity invariance for 
quark spin orientations perpendicular to the q-C plane and is allowed by time-reversal 
invariance due to the soft interactions of the fragmenting quark with external fields. 

This idea was exploited in Ref. 5 to explain the single spin asymmetries observed in 
p^p — i> ttX processes i and where, essentially, it is assumed that parton c is produced 
in the forward direction and the final hadron pt is due to its transverse k± inside the 
jet. One cannot expect such a model to work at large vt., as the results clearly show. 

Another possible k± effect, suggested by Sivers Q'B, may originate in the dis- 
tribution functions. Similarly to the Collins effect in the fragmentation process a 
dependence on the hadron spin may remain in the k± dependent quark distribution 
function fg/At{x, k±), so that the difference 

A^/g/AT {X, k^) = /g/^T (X, fcx) - fg/Al {x, k^) 

= fg/AT{x,k^) - fg/Al{x,-k^) (10) 
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can be different from zero. fq/A^ {x, k±) is the number density of partons a witli mo- 
mentum fraction Xa and intrinsic transverse momentum fe^ in a transversely polarized 
hadron A. 

Taking intrinsic k± into account Eq. (|^) becomes 



c 



E 

a,b,c,d 



d'^k I dx„ dxh 



TIZ 



da 



X fa/A^{.Xa, k±) fb/sixb) "^(^i) Dc/ciz) 



so that, adopting obvious shorter notations 



d'^k±dxa dxb 



TTZ 



(11) 



(12) 



fa/A^iXa, k±) - fa/A^iXa, -fc±)J fb/B{Xb) da{k±) Dc/c{z) • 

The above new quantity (0), divided by twice the unpolarized fc^ dependent 
distribution function 2fq/A{xa, k^) = fq/A^{xa,k^) + fq/A^{xa,k±), can be regarded 
as a single spin asymmetry or analysing power for the ^ a + X process; similarly 
to the quark fragmentation analysing power, Eq. (P), it may be different from zero 
for hadron spin orientations perpendicular to the A-q plane and taking into account 
initial state interactions between the two colliding hadrons. This quantity plays, for 
single spin asymmetries, the same role plaid by distribution functions in unpolar- 
ized processes. If we define the polarized number densities in terms of distribution 
amplitudes as 



fa,XJA^ {Xa,k±a) 

then we have, in the helicity basis 

i/AT 



I^Ax,,A„;T(^a' k±a)\'^ 



(13) 



^^fa/A^ iXa,kj 



E [Gx^^,X^,+ {Xa,k^a) Gx^^X;4Xa,k^a] 



(14) 



which shows the non diagonal nature, in the helicity indices, of /^^(xa, k±a)- 
Notice that, due to the odd k± dependence of /+(^, one has 



Jd'k^rJ^{k^)^{k 



d'kirj^ik^ 



{kj_)^>o 



da 
dt 



da 
di 



-k 



(15) 
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where da / di(k±) — da / di{—k±) = 0{k^/pT)i which clearly shows that this is a higher 
twist effect. 

A simple phenomenological model was developed in Ref. 9 for the p^p — ^ vrX 
process by assuming that the dependence of on k±a is sharply peaked around an 
average value /c°a- 

lt{Xa.k^) = (l^-l - A^aX-(l - X^f^ (16) 

with, from a fit to the data 0, 

^ = 0.47a;°■^«(l-x)°•^^ (17) 

where M is a hadronic mass scale, M ~ 1 GeV/c. 

Inserting Eqs. ([T^), ([T5|) and (p^ into Eq. (|D yields 



^ _ 1 EA^a/^a^aC^Xfcfc°,x°°(l-xJ^"/,/p(x,)[rfa(A:U-^'»(-feU]-Pvc(^^ 
^ M EJdXadxb fa/piXa) fb/p{xb)daD^/c{z)/z 

(18) 

In order to give numerical estimates of the asymmetry (|18D we have taken fq^q^g/p 
from Ref. 11, D^^jq^q from Ref. 12 and D^^jg from Ref. 13. Given the very limited 
Pt range of the data we have neglected the QCD Q"^ dependence of the distribution 
and fragmentation functions. We have only considered contributions from u and d 
quarks inside the polarized proton, which certainly dominate at large xp values, that 
is a = M, (i in the numerator of Eq. ([181) . Instead, we have considered all possible 
constituents in the unpolarized protons, with k±_ = 0, and all possible constituent 
fragmentation functions. 

By using Eq. ( p!7D into Eq. ([l^), together with the unpolarized fa,b/p and D^^/c 
functions, we remain with an expression of still dependent on a set of 6 free 
parameters, namely iV^, and Pa (« = u,d), defined in Eq. (0). We have obtained 
a best fit to the data ^, shown in Fig. 1, with the following values of the parameters: 

Na eta Pa 

u 5.19 2.10 3.67 , , 

(19) 

d -2.29 1.43 4.22 



As the experimental data a cover a px range between 0.7 and 2.0 GeV/c we have 
computed Ajq at a fixed value px = 1.5 GeV/c. 

Notice that the above values (|19]) are very reasonable indeed; actually, apart 
from an overall normalization constant, they might even have been approximately 
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Fig. 1. Fit of the data on Ajv i, with the parameters given in Eq. ([Tg|); the upper, 
middle and lower sets of data and curves refer respectively to tt^ , tt*^ and tt~ . 



guessed. The exponents au,d and f3u,d are not far from the very naive values one can 



obtain by assuming, as somehow suggested by Eqs. (|T^) and ([T^, that /+'^(xa 



-/p,+{^a), where fa,+{-)/p,+ {^a) denotes, as usual, the number density 
of quarks with the same (opposite) helicity as the parent proton. Also the relative sign 
and strength of the normalization constants A^^ and A*",;; turn out not to be surprising 
if one assumes that there might be a correlation between the number of quarks at a 
fixed value of k±a and their polarization: remember that, according to SU{6), inside 
a proton polarized along the y direction, Py = 1, one has for valence quarks Py = 2/3 
and P;j = -1/3. 

It might appear surprising to have approximately opposite values for the vr"*" and 
7r~ asymmetries, as the data indicate, and a large positive value for the 7r°; one might 
rather expect Ajv ~ for a 7r°. However, this can easily be understood from Eq. ([T8|) 
which we simply rewrite, for a pion vr*, as = A/"*/©*, if one remembers that from 
isospin symmetry one has: 



(20) 



Eq. (|T^) shows that the relation ( pO| ) also holds for Af^ and so that 



4° 



AT^+AT 
V+ + V- 



1 + 



1 + 



(21) 
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It is then clear that A]^ ~ iniphes ~ only if ~ i.e. if the 

unpolarized cross-sections for the production of a 7c~ and a tt^ are approximately 
equal. This is true only at ~ 0. At large xp the minimum value of Xa kinematically 
allowed increases and the dominant contribution to the production of tt"*" and vr" 
comes respectively from fu/pixa) and fd/p{xa) [see the denominator of Eq. 
It is known that fd/pi^a) / fu/p{xa) when Xa — ^ 1; this implies that T>^ /T>^ 
decreases with increasing xp-, so that at large xp we have jD^ ^ 1 and ~ 
A'^{\- 2V-/V+) ^ A%. Such a trend emerges both from the experimental data and 
our computations. 

We have discussed how to compute single spin asymmetries in large pt inclu- 
sive production within the framework of the factorization theorem and perturbative 
QCD, showing that single spin effects of order k^/px can be related to non pertur- 
bative intrinsic properties of quark fragmentations and distributions, i.e. the quark 
fragmentation analysing power and the quark distribution analysing power, the former 
was first suggested by Collins i and the latter by Sivers 011. We have shown a possible 
description of the single spin asymmetries observed in p^p —>■ nX via the intrinsic k± 
effects in the quark distribution functions. 

A definite test and a better evaluation of these non perturbative properties requires 
further and more refined applications of the same idea and more theoretical work. In 
particular one might consider the following processes: 

- p^p — i> ttX; a straightforward application of the model previously described, with 

the same set of parameters, gives results comparable with some existing data 

E. 

- p^p — i> 7X; in such there cannot be any fragmentation effect and one 

should be able to single out the importance of the Sivers effect and the quark 
distribution analysing power. 

- ip' — > hX; in such the inclusive production of a hadron in DIS, there 

should be no effect from the quark distribution analysing power, as any initial 
state interaction would be of higher electromagnetic order and negligible. It 
should then allow an evaluation of the Collins effect. 

- p^p vrX, with a correct inclusion of both the Collins and Sivers effect. 

Let us finally mention that we have assumed the validity of the QCD factorization 
theorem also at higher twist. This has been discussed in the literature: an approach 
similar to that discussed here i has been advocated, in the operator language, by 
Qiu and Sterman who use generalized factorization theorems valid at higher twist 
and relate non zero single spin asymmetries in pp collisions to the expectation value 
of a higher twist operator, a twist-3 parton distribution, which explicitly involves 
correlations between the two protons and combines quark fields with a gluonic field 
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strength. However, they still consider only coUinear partonic configurations so that, 
in order to obtain non zero results, they have to take into account the contributions 
of higher order elementary interactions. Some more theoretical work is still necessary. 
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